A Parallel Three-Dimensional Incompressible Navier-Stokes Solver with a Parallel Multigrid Kernel

The development and applications of a parallel, time-dependent, three-dimensional incompressible Navier-Stokes flow solver and a parallel multigrid elliptic kernel are described. The flow solver is based on a second-order projection method applied to a staggered finite-difference grid. The multigrid algorithms implemented in the parallel elliptic kernel, which is used by the flow solver, are V-cycle and fidl V-cycle schemes. A grid-partition strategy is used in the parallel implementations of both the flow solver and the multigrid kernel on all fine and coarse grids. Numerical experiments and parallel performance measurements show the parallel solver package is numerically stable, physicrdly robust and computationally efficient. Both the multigrid kernel and the flow solver scale well to a large number of processors on Intel Paragon and Cray T3D/13E for two and three dimensional problems with moderate granularity. The solver package has been carefully designed and implemented so that it can be easily adapted to solving a variety of interesting scientific and engineering flow problems. The mde is portable to parallel computers that support MPI, PVM and NX for interprocessor communications.